library(tidyverse)
library(DeclareDesign)
library(rdss)
library(CausalQueries)

sims <- 2000
bootstrap_sims <- 2000

source("code/declarations/declaration_16.1a.R")
source("code/declarations/declaration_16.1b.R")

diagnosis_16.1 <-
  diagnose_design(
    declaration_16.1,
    sims = sims,
    bootstrap_sims = sims,
    make_groups = vars(XY)
  )

write_rds(diagnosis_16.1, file = "diagnosis_objects/diagnosis_16.1.rds")

